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1 Introduction 

Hendrik Brugt Gerhard Casimir (1909-2000) was born in The Hague (the Nether- 
land). He studied physics and obtained his Ph.D in 1931 at the University of Leiden. 
He visited Bohr in Copenhagen and was an assistant to Pauli at Zurich in 1938. Then 
he became a physics professor at Leiden University. In 1942, during world war II, 
he moved to the Philips Research Laboratories in Eindhoven where he remained 
an active scientist. He became a member of the board of directors of the Philips 
company in 1956. In addition to his professional life in industry, Casimir was a 
deep and broad-minded theoretical physicist. He was active in pure mathematics 
(Lie groups), in low temperature physics, irreversible phenomena, and fluctuation 
induced forces (the so-called Casimir effect) 

The seminal papers that lead to the Casimir effect where motivated by discrep- 
ancies between theory and experiments on colloidal particles in suspension. In 1947 
Verwey and Overbeek, applying the London-van der Waals theory to their obser- 
vations on the attraction between large particles in suspension, concluded that this 
theory was not in agreement with their experimental results. The discrepancy could 
be removed if the potential between the particles decreased more rapidly than r~ 6 . 
They suggested that the influence of retardation on the interaction should be taken 
into account. In [1], Casimir and Polder indeed found that the finiteness of the speed 
of light causes the asymptotic potential to decay as r~ 7 . The (ultimately simple) 
result follows from an elaborated 4th order perturbation calculation in quantum 
electrodynamics and will not be presented in these notes. Looking for a simpler 
derivation, Casimir discovered [2] that the change of the zero point quantum energy 
due to the presence of two metallic plates generates a macroscopic force between the 
plates. That was a striking macroscopically observable manifestation of the effects 
of vacuum fluctuations in quantum electodynamics. Then Casimir applied similar 
ideas to the derivation of the retarded van der Waals forces [3] that we will discuss 
in section 4.3. 



2 



Nowadays the field of Casimir forces has grown immensely wide. The Casimir 
effect could be tentatively defined as the change of the vacuum energy (zero point 
energy) by external constraints and its consequences. In condensed matter physics it 
covers the field of induced forces by vacuum and/or thermal fluctuations on atoms, 
molecules and macroscopic bodies. One can distinguish 

• Electromagnetic Casimir effect 

at the macroscopic level, the forces between conductors and dielectric bodies; 
at the atomic level, dispersion forces, van der Waals forces, molecular attrac- 
tions, atom-walls interactions, 

• Casimir effect in general quantum field theory 

modification of the vacuum energy in presence of external sources or geomet- 
rical constraints, 

• Casimir effect in particle physics 
e.g. in the bag model of hadrons, 

• Casimir effect in cosmology 

e.g. modification of the zero point energy in topologically non trivial spaces 
compared to Minkowski space, 

• Casimir effect in critical phenomena 

forces between boundaries and layers due to long range order at critical points 
and in phases with broken continuous symmetry, 

• Dynamical Casimir effect 

time dependent Casimir forces generated by moving boundaries or constraints. 

The large number of text books and extended reviews (not mentionning research 
papers) bears testimony for the activity and broadness in the subject. An non 
exhaustive list, focusing on one or the other aspect of Casimir forces, can be found 
at the end of these notes [4], [5], [6], [7], [8], [9], [10], [11], [12], [13], [14]. 

It is obviously not possible to cover the whole domain in six one hour lectures, and 
we shall therefore mainly concentrate on the traditionnal electromagnetic Casimir 
effect. We will however focus attention on a more detailled and complete statistical 
mechanical treatment than that usually found in the above texts, hoping to show 
that some new information can still be brought to this by now old and venerable 
subject. 

In sections 2.1-2.2 we review the standard treatment of the Casimir force be- 
tween two metallic plates at zero temperature. The plates are caracterized by the 
boundary conditions of macroscopic Maxwell fields at metallic interfaces, namely 
the vanishing of the part of the electric field tangent to the surface. This geomet- 
rical constraint modfies the field eigenmodes depending on the distance d between 



3 



the plates. The d-dependance of the modified zero point energy is the source of the 
Casimir force. The generalization of Casimir's calculation to non zero temperature 
is presented in section 2.3. Each field mode is now a thermalized quantum mechan- 
ical oscillator with frequencies obtained from the previously described macroscopic 
boundary conditons. The total force between the plates, due to the sum of vacuum 
and thermal fluctuations of the field, is normalized in such a way that a single plate 
immersed into an infinitely extended radiation field does not experience any force. 
It is of interest to note that the high temperature asymptotics of the force is purely 
classical i.e. independent of the Planck constant. This is the classical regime of 
the Casimir force. We stress again that all these considerations are based on the 
premises that the conductors are treated as macroscopic bodies without internal 
fluctuations, namely the microscopic charges and field fluctuations inside the con- 
ductors are ignored. Only the photons outside the conductors are subjected to a 
statistical distribution. In this situation we refer to conductors as being inert or 
dead. 

In sections 3.1-3.2, we study a purely classical model for the Casimir force [15]. 
Classical charges are confined into two globally neutral slabs separated by a distance 
d. All the charges are in thermal equilibrium and interact by means of the pair- 
wise Coulomb interaction. We then proceed to a direct calculation of the average 
Coulomb force by unit surface using the standard methods of classical statistical 
mechanics of fluids (integral equations and Mayer graph summations). We find that 
the asymptotic force is half the value of that found following Casimir's method based 
on dead conductors. This confirms the result of alternative calculations previously 
made for similar classical models [16], [17]. In these works the premises are very 
different : the field fluctuations (here of purely electrostatic nature) are generated 
by those of the charges in the conductors. When the microscopic charge fluctuations 
inside the bodies are taken into account we will speak of living conductors. 

The factor 1/2 of discrepancy between the Casimir force for dead and living con- 
ductors in the classical regime could be attributed to the fact that in purely classical 
models, the transverse part of the electromagnetic field responsible for the radiation 
does not come in. This leads us to to consider a more complete model where the force 
between two quantum plasma slabs is computed in the framework of non relativistic 
quantum electrodynamics including both quantum and thermal fluctuations (section 
3.3). The analysis, based on the path integral representation of the quantum gas 
[18], shows that in the semi-classical regime, quantum corrections do not alter the 
asymptotic form of the force previously found in the classical models. The quantum 
effects manifest themselves as tiny corrections of diamagnetic type occuring at the 
subdominant order. The conclusion (see a more thorough discussion in section 3.3.3 
and in [19], [20]) is that the treatment of conductors as dead, prevailing in the whole 
literature, is not physically correct as soon as the temperature is different from zero. 
Fields and particles do fluctuate in the conductors which make the enforcement of 
inert boundary conditions inadequate. 

The section 4 is devoted to some aspects of dispersion forces. After briefly re- 
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calling the standard non retarded van der Waals-London forces between two atoms 
in vacuum, we report in section 4.2 on the status of these forces in a dilute medium 
having non zero temperature and density [21]. Here the situation is not quite ele- 
mentary because one has to treat in a coherent way both the quantum mechanical 
binding leading to the formation of atoms and the collective screening effects that 
are always present as soon as there is a fraction of ionized charges. In the proper 
scaling regime called the atomic limit (high dilution and low temperature) we can 
give the exact asymptotic form of the correlations up to exponentially small terms 
as T — ■> 0. One finds that large distance atomic correlations reduce to the standard 
van der Waals-London potential r~ 6 with some polynomial corrections in tempera- 
ture inherited from collective screening effects. It turns out that unbound charges 
and atom-charges both undergo van der Waals type of interaction as a consequence 
of the fact that correlations between quantum charges always have a r~ 6 decay 
(exponential screening never holds in quantum plasmas). 

Following Casimir's orginal ideas about the effects of vacuum fluctuations, we 
present in section 4.3 an elegant derivation of the retarded van der Waals forces 
at zero temperature. The derivation relies on the assumption that in addition to 
the electromagnetic field radiated by the sources there is always an underlying free 
quantum electromagnetic field present whose vacuum fluctuations is at the origin of 
the van der Waals forces. 

Finally section 4.4 gives a short account of forces between dielectric bodies es- 
sentially following the famous Lifshitz semi-macroscopic theory. Details and devel- 
opments can be found in text books. 

In section 5 we touch the subject of the Casimir effect in critical phenomena by 
considering the simplest quantum system that exhibits a phase transition, the free 
Bose gas [22]. It is shown that the grand-canonical potential of a Bose gas in a 
slab at the critical value of the chemical potential has finite size corrections of the 
standard Casimir type. They can be attributed to the existence of long range order 
generated by gapless excitations in the phase with broken continuous symmetry. 

2 The historical calculation 

2.1 The metallic cavity 

Here we follow partly [4], sec. 2.7, and [5] (B. Duplantier), see also [14]. 

We consider an empty cubic box A with perfectly conducting boundaries. It has 
thickness d and lateral sizes of surface L 2 : A = {x = (x, y, z) |0 < x < d, — -§<?/< 

— -§ < z < -§}. The electric field in the box is solution of the Maxwell equations 
in empty space (Gauss units) 

V 2 E(x,t)- 1^1^=0, V-E(x,t)=0 (1) 
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with the boundary conditions 

E tg (x,i)=0, xe<9A (2) 

since the tangential components of the field E tg have to vanish on perfectly conduct- 
ing walls. The eigenmodes are found by setting E(x, t) = Re(E(x, uj)e~' tuJt ) in (1) 
leading to the eigenvalue equation for the complex amplitudes E(x, a;) (the Helmotz 
equation) 

V 2 E(x, u) = -^ E ( x - v • E ( x > ^) = °> E tg( x > = 0, xe9A. (3) 

The solutions, labelled by wave numbers k and polarization indices A, are 

E k , x e x (\) cos(k x x) sin(k y y) sin(k z z) 
Eic,a(x) = I E KX e y (X)^^sm(k x x) cos(k y y) sin(k z z) (4) 
E KX e z (X)^J^sm(k x x) sm(k y y) cos(k z z) 

where ->/8/|A|, |A| = L 2 d, is a normalization factor, e(A) {e(A) -k = 0, e(A) -e(A') = 
^a,a')A = 1,2)} are two polarization unit vectors othogonal to the wave number k, 
and E^x is a complex amplitude. The wave numbers and the eigenfrequencies are 
of the form 

7 _ T^x_ , _ ^hs_ h _ ™^ - n 1 9 

K X j Ky [j 1 Xj ' ^ Z ' ' ' ' ' ' 



c<j k = c|k 



The corresponding magnetic field amplitude B(x, a;) is obtained from Faraday equa- 
tion B(x, uj) = — i^V AE(x, uj). One sees on (4) that if one of the integrers n x ,n y , n z 
vanishes, E kiA (x) is necessarily directed in a single direction, i.e. there is only one 
possible polarization state in this case. The energy of one mode is 



i- [|Re(E k , A (x))| 2 + |Re(B k , A (x))| 2 ] 

1 / 77 1 * 77" I 77" 77 1 * s 

/ -^k^-^k^ + -&k,A-fr k ,A 



8vr 



(6) 



The electric and magnetic energy contribute the same amount and we have kept the 
order of the products as they occur in the calculation. 

The quantized electric field (still noted E kiA (x)) is obtained as usual by replacing 
the classical amplitude 

E KX -> (y^Khu^j a M , ££ >A -> [y^KhJ'kj a k ,A 



by the (dimensionless) annihilation and creation operators of photons with the com- 
mutation relations [ak,\,at' \>] = ^k,k'^A,A'- 

Then the quantized energy of one mode takes the form 

A,k,A = ^k ( 2 ) = Mc |a kA a k , A + - I (7) 

Finally the total energy results from the summation over all modes k, A is (taking 
into account that the mode eigenfunction in (4) are orthogonal and normalized) 

H A = H AXX = hw k a kiA a k ,A + ^ • (8) 

k,A k,A k,A 

The notation means that only one polarization is possible when one of the wave 
numbers k x , k y , k z is equal to zero. The (infinite) last term of (8) represents the zero 
point energy of the quantum electromagnetic field. There are two ways of thought 
about it. It merely appears in (8) as an additive constant to the total photon energy 
so that one can dispense with it by redefining the energy of the photon vacuum 
state |0 > to be equal to zero. This is commonly done by writing the creation and 
annihilation operators in (7) in normal order. More fruitful is Casimir's view that 
the vacuum energy 




also represents mean square fluctuations of the fields in the box A that exist even 
in the absence of photons. These pure vacuum fluctuations may have physically 
observable effects since they depend on the geometry (shape, size) of the spatial 
domain that constrains the field. 

2.2 Force between macroscopic metallic plates induced by 
vacuum fluctuations 

We are interested in the force by unit surface f(d) = F(d)/L 2 induced by vacuum 
fluctuations between two faces of the metallic box at distance d. Since there are only 
two fondamental constants c and h available in the problem, dimensional analysis 
shows that f(d) ~ hc/d 4 , the point being to determine a finite proportionality 
coefficient. There are several ways to regularize the infinite energy (9) to extract 
physically meaningful quantities. A perfect metal has a static dielectric constant 
e = oo. One argument consists in observing that a real metal is characterized by a 
frequency dependent dielectric function e(u) such that e(uS) — > oo, uj — > 0, but which 
tends to the vacuum value eo as uj — > oo, namely when uj ^> uj a , uj a a characterisitic 
atomic frequency. Hence high frequencies should not contribute to the force, and 



7 



for this reason one introduces a cut-off function g(uj/u a ) in (9) such that g(0) = 1 
and g{u /u) a ) — > as uj — > oo 



k,A 



(10) 



The cut-off function will be removed at the end of the calculation by letting uj a — > oo. 

Extending the plates to infinity in the y, z directions gives the energy per unit 
surface 



w vac (d) 



lim -| 



7T Z 



1—\ dk y dk z - V + A; 2 

^0 ./o ^„ =0 vv 



:n) 



where q is the two dimensional wave number vector in the (y, z) plane, q = |q| and 
the second line results from the integration in polar coordinates in this plane with 
angular sector 2n/4. The prefactor 2 is due to the two polarization states and 
means that the term n = must have an additional factor 1/2. With the change of 
variable q — > uj = u> n (q), ujduj = c 2 qdq, u vac (d) can also be written as 



u vac (d) 



1 



2ttc 2 

hen' 2 
~2~¥ 



E 

n=0 



poo 


hu g ( — ) 


J dujuj 


' cirn/d 





5> 



r?. 



n=0 



where we have defined 



dvv 2 g 



TIC 
dWn 



(12) 
(13) 

(14) 



Finally the work required to bring the plates from a large separation d ~ oo to the 
actual separation d and the force between the plates are 

w (d) = u vac (d) - lim w vac (rf), r c {d) = -77iu vac (d) + lim -^zu^d) . (15) 

d^oo (7(2 a!^oo (7(i 



One sees on (13) that 
u vac (d) ~ 



1 d 



27CC 3 7T 



2d 3 7o 



rfsF(s) . 



hu g [ — 



OJr, 



, d — > oo 



(16) 
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Using the above definitions w(d) can be cast in the form 



w(d) 



7T 2 He 

~2dF 



00 1 poo 

J2 F (n) + -F(0)- dsF(s) 

n=l J ° 



(17) 



In (17) we have explicitly singled out the n = term of the Y^' n sum. The asymptotic 
evaluation of the difference between the sum and the integral (17) is provided by an 
application of the Euler-MacLaurin theorem [23], [24] Sec. 6.6.4 



00 poo 1 r> 

5>(n)= / dkF(k) - -F(0) - -f F« 

n=l J° 



(0) 



f^)(0)-f^)(0) + ... 



(18) 



with 5j the Bernouilli numbers B 2 — |, -64 = 
easily calculated from (14) 



^, The derivatives F^(s) are 



FW(s) 
F (1) (0) 



7TC 



0, F (3) (0) = -2 + 



dw a 



F (j) (0) = C 



dw a 



, 3 > 5 ■ 



(19) 



This holds under the assumption that the cut-off function verifies g(0) = 1, all its 
derivative are finite at the origin and vanishing at infinity. With (18) and (19) and 
letting the cut-off ui a — > 00 we obtain Casimir's result: 



w(d) 



7T 2 hc 

720rf 3 ' 



f™(d) 



7T 2 hc 

"240rf 4 



(20) 



In this calculation, the electromagnetic field is entirely enclosed in the cavity, the 
outside of it being void of matter and of electromagnetic energy. We can take the 
alternative equivalent view that there is also a field in the space external to the cavity. 
Then the external face of the plate at d will be subjected to a force in the opposite 
direction due to the vacuum fluctuations in the semi- infinite space to its right, 
namely /™ c = - (-lim^oo ^w vac (d)), so that the total force —§- d u vac {d) + f™(d) 
is the same as (15). This amounts to normalize the force so that a single plate in 
infinite space feels no resulting force. 



2.3 Force induced by the thermal fluctuations 

If the walls of the cavity have a certain temperature T ^ 0, the photons inside the 
cavity are thermalized at the same temperature by their interactions with the atomic 
matter constituting the metal. Assuming that the thermalization process has taken 
place, the photon-atoms interactions are no more described at the microscopic level 
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in the standard treatment of blackbody radiation, but replaced by the macroscopic 
boundary condition (2) of Maxwell fields at a metallic interface. Hence, each mode 
of the field behaves as a quantum mechanical oscillator at temperature T. 

When temperature is introduced in this way there is a new typical length in the 
problem, the thermal wave length of the photon /3hc providing the dimensionless 
parameter 

Birhc . . 

" = — ■ < 21 > 

The large values of a correspond to low temperature or short distances where the 
quantum aspects of the electromagnetic field play a dominant role. The small values 
of a (high temperature or large distance) will yield the classical limit. 

The free energy of one photon of frequency with energy levels Hcu^n, n = 
0, 1, 2, . . . is (disregarding the zero point energy huk/2 at the moment) 

-/T 1 In (jr e ~^A = /r 1 *^), X (v) = In (1 - e-") . (22) 

\n=0 / 

Hence the thermal free energy of the electromagnetic field in the metallic box is 

$t,a = Y^P^xiPfu^) . (23) 

k,A 

The free energy per unit surface between two faces at distance d is obtained by 
performing exactly the same steps as in (11)— (13) replacing there ^fkug(uj/u a ) by 

rr OO /*00 

= -jt = wE' duJUJ WM] • 

(24) 

The radiation pressure due to the thermal photons inside the cavity on the plate at 
d can be written in terms of the parameter a as 

#V) = -TgpM = E(-) 2 X'(-) ■ (^) 

n=l 

This pressure is positive since xi v ) is negative. 

The pressure on the plate due to thermal radiation in an infinite half space is 
given by the well known formula 

p-looHJunp-M^^. (26) 

It is obtained by taking the continuum limit of the sum (25) as a — > 

pS^oo) = - Jim — - 2 / dvv 2 In (l - e~ v ) (27) 
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The result (26) follows when we insert the expression of a (21) and use the relations 

poo 

I p = dvv p ln(l-e- v ) =p\ CGo + 2) (28) 
Jo 

with ( the Riemann function 

CM = £^ «<) = =;. (29) 



nP * w 90 

n=l 



2.3.1 Short distance or low temperature limit 

The large a expansion is obtained by noting that x( na ) ~ ~~ e~ na , a — > oo. Keeping 
the dominant term n = 1, (25) gives 

^ d(rf) = rf^ ^ + • (30) 

Thus the pressure of a very thin black body is exponentially small at fixed T (or at 
fixed d in the low temperature limit). 

2.3.2 Long distance or high temperature limit 

In order to obtain the small a expansion of the sum in (25) one uses the Poisson 
summation formula which in our case reads 1 2 

_. oo oo 

-F(0) + ^F(n) = 7rC(0) + 27r^C(27rn) (31) 
2 

n=l ra=l 

where C(k) are the Fourier coefficients of F(s) 

1 f°° 

C(k) = - / dsF(s) cos ks . (32) 
71" Jo 

With F(s) = s 2 x{as) = s 2 ln (1 - e~ as ) one finds 



1 d 

C(k)- 



IT 



dk 2 



k 2 k \ « 



(33) 



The value of C(0) = \ / °° dss 2 In (1 — e _QS ) = — 4^ can be found from (28) or 
simply by expanding cothx for small argument. For k = 2nn, n 7^ one uses 
cothx = 1 + 0(e~ 2x ) leading to 

C(2irn) = —1— - ,l a N + O fexp (-b-)) , b>0. (34) 
v ; (27m) 3 7r(27rn) 4 V F V a) J y ' 



1 This is the version of Poisson formula adapted to the evaluation of sums on positive integers. 
2 The Euler-MacLaurin formula is not adapted here because the derivatives of the function 
In (1 — e~ v ) diverge at v — and therefore an exact summation is needed. 
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This is inserted in Poisson formula (31) with F(0) = and the definition (29) of the 
Riemann function 



and finally 

^-w-S+s£+°M-!)) < 36 > 

This is the radiation pressure (26) together with additional high temperature or finite 
distance corrections. One notes that the d~ A term has a temperature independent 
coefficient which is exactly the same as that of the vacuum Casimir force (20) in 
magnitude but with the opposite sign. 

2.3.3 The total force 

The quantity of interest is the total force f(d) on the plate at d when it is also 
submitted to black body radiation from the infinite half space to its right producing 
the pressure —p^ d (oo). Thus adding also the vacuum Casimir force that was not 
taken into account in the above thermal calculations one has 

f(d) = p? d (d) - p? d (oo) + r c (d) . (37) 

This amounts to consider that a single plate immersed in an infinitely extended 
radiation field experiences no net force and f(d) is the force that develops when a 
second plate is present at finite distance d. From (20), (26), (30) and (36) one finds 
the large and small a expansions of f(d) 

K 2 hc 7T 2 1 7T „ , „ N . . 

'w=-^-«^W (e h a ~°° (38) 

and 

W-Me +0 (^H))- a - ' (39) 

Concerning the latter result there is a number of remarkable points: the dominant 
order d~ 3 has a purely classical expression independent of Planck's constant and it 
is still attractive; moreover at next order d~ 4 there is an exact compensation of the 
Casimir vacuum force. 

2.4 Some experimental results 

We first discuss the order of magnitude of the Casimir force. For a separation 
d ~ 10 _6 m = 1/im and at room temperature T = 300-fT one has a = /3nh~c/d ~ 24. 
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We are thus in the large a regime and the temperature corrections (last term of the 
r.h.s of (38)) have a factor e~ 24 !. Then the ratio of the blackbody pressure (second 
term of the r.h.s of (38 ) to the Casimir force (20) is lQn 4 a 4 /3 ~ 1.5 1CT 3 so that 
vacuum fluctuations dominate the black body pressure. One can therefore observe 
a macroscopic electromagnetic force between the plates at room temperature entirely 
du to vacuum zero point energy ! 

Direct experiments in the planar geometry are difficult because of a number of 
perturbations and have only be performed recently with sufficient accuracy. The 
parallelism of the plates has to be perfect, corrections due to to additional elec- 
trostatic charges, to the roughness of the surfaces and to their finite conductivity 
have to be taken into account. We quote the recent results of G. Bressi et al. [25]. 
The agreement with the d~ 4 decay is excellent. Moreover the calculated value of 
the Casimir amplitude K c = n 2 hc/2A0 is K c = 1.3 1CT 27 Newton m~ 2 whereas 
these authors have measured the value Kc = 1.22 ± 0.18 10~ 27 Newton m~ 2 . One 
concludes that the Casimir formula (20) for the zero temperature force is now ex- 
perimentally validated. 



2.5 Non planar geometries 

The calculation of the Casimir force for other types of geometries is not easy because 
of the need to determine the electromagnetic eigen modes with metallic boundary 
conditions on non planar surfaces. Earlier experiments where performed on the 
force between a sphere and a plate. With this arrangement one avoids the problem 
of contolling the parallelism of plates. A common approximation for the force is 
provided by the Derjaguin construction [4], [6]. It amounts to decompose the sphere 
into a succession of concentric cylindrical shells with axis perpendicular to the plate. 
The force between each of the cylinder is assimilated to the Casimir planar formula 
(20), and the total force on the sphere is recovered be summation of all these force 
elements. Within this approximation the zero temperature force exerted on the 
sphere by the plate is 

^/ , n 3 he R , . 

F ^ R ^-lm^ (40) 

where d is the distance between the surface of the sphere and the plate, R is the ra- 
dius of the sphere and i? S> d. This formula has received a rather good experimental 
verification [26]. 

An important issue is the sign of the Casimir force. One has seen that the force 
between two plates or a plate and a sphere is attractive. Motivated by the Casimir 
model of the electron, the question of the vacuum energy force on a conducting 
spherical shell was raised. In this model the electron is regarded as a uniformly 
charged conducting shell of radius r with total charge e. The electrostatic repulsive 
energy of the shell e 2 /2r could be balanced by an attractive Casimir energy which 
should be of the form of the form — C hcjlr for dimensional reasons. This con- 
figuration would be stable if the dimensionless constant C is positive and has the 
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exact value C = e 2 /he. With this beautiful idea the electron constitution would be 
entirely explained on the basis of electromagnetic effects. After several works, it was 
eventually established by Boyer [27] that in the situation of a spherical shell, the 
Casimir force is repulsive (it tends to dilate the sphere), so invalidating Casimir's 
proposal. 

Balian and Duplantier have treated the problem of a number of perfectly con- 
ducting shells of general shapes [28] . The shells are idealized as smooth surfaces on 
which the field has to satisfy the macroscopic metallic boundary conditions. The 
authors provide the general form of the free energy at all temperatures and depend- 
ing on the topology of the surfaces. The eigenmodes enter through the density of 
states, itself given by the spectral function associated to the electromagnetic Green 
functions. The latter can be calculated in a perturbative scheme called "multiple 
scattering expansion". Previously known particular cases corresponding to simple 
geometries are recovered, and results for more general shapes are obtained in terms 
of curvatures and genders of the surfaces. A description of this mathematically nice 
and general theory can be found in [5] and [14]. One should remember that all the 
results quoted in this section were based on the hypothesis of "inert" conductors, 
i.e. neglecting charge and field fluctuations inside the conductors. We come to this 
point in the next section. 

3 Statistical theory of the classical Casimir effect 
3.1 The model 

Our goal is to analyse the Casimir force at a more fundamental level by introducing 
in the description the microscopic degrees of freedom of the atoms constituting the 
plates. For simplicity, and motivated by the high temperature result (39), we shall 
first present a purely classical model. The transverse part of the electromagnetic 
field is not considered here and the forces are purely electrostatic. This model is 
treated in detail in [15]. 

The two plates A and B consist of classical point charges confined to two planar 
slabs Aa(L,o) and Ab(L, b) in three-dimensional space. One can think of them as 
weakly coupled plasmas or electrolytes that are conducting in the sense that they 
are characterized by a microscopic screening length £ D which is of the order of the 
interparticle separation. 

The slabs have thickness a and b, surface L 2 , and are separated by a distance d : 

A A (L, a) := {r = (x,y) \xE [-a, 0], y G [-§ , § ] 2 } 

A B (L, b) := {r = (x, y) | x G [d, d + 6], y G [-f , § ] 2 } (41) 

and we shall assume that all the lengths a, b, d are much larger than Id. 

Plasma A (B) is made of charges e a (ep) of species a G Sa (P G Sb) where 
Sa and Sb are index sets for the species in Aa(L,q) and As(L,b) respectively. We 
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assume both plasmas to be globally neutral, i.e. carrying no net charge, 

a b 

where J2 a (J2b) extends on all particles in Aa(L, a) (A B (L, b)). For a particle located 
at r we will use the generic notation (7 r) where 7 G Sa if r G Aa(L, a) and 7 G 
if r G A#(L,ft). The space external to the slabs is supposed to have no electrical 
properties, its dielectric constant being taken equal to that of vacuum. The charges 
are confined in the slabs by hard walls that merely limit the available configuration 
space to the regions (41). 

All particles interact via the two-body potential 

^(7,7 / ,|r-r'|) =e 7 eyu(r-r / )+u SR (7,7 / ,|r-r / |), (43) 

where v(r — r') = l/|r — r'| is the Coulomb potential (in Gaussian units) and 
^sr(7, 7') r— r') is a short-range repulsive potential to prevent the collapse of opposite 
charges and guarantee the thermodynamic stability of the system. 

The total potential energy U consists in the sum of all pairwise interactions, 
separated into three contributions according to whether they take place between 
two particles of A, of B, or between a particle of A and a particle of B : 

U = U A + U B + Uab- (44) 

On the microscopic level, the force between configurations of charges in the two 
plasmas is the sum of all pairwise forces exerted by the particles of A on the particles 
of B : 

F Aa _ As := - 

a b 

r a eA A (L,a), r b eA B (L,b) (45) 

and F SR is the force associated to the short-range potential i>sr- For simplicity we 
assume that the range of i>sr is finite so that FsR(o;a, A>, r a — r&) vanishes as soon 
as d is large enough, and we will omit it in the following. 

Both plasmas are supposed to be in thermal equilibrium at the same temper- 
ature T. The statistical average (•••).£ is defined in terms of the Gibbs weight 
exp(— (3U), (3 = (ksT)' 1 , associated with the total energy (44). There is no need 
to explicitly specify the ensemble used here (canonical or grand canonical) provided 
that the global neutrality constraint (42) is taken into account. The average parti- 
cle densities p^^r) are expressed as averages of the microscopic particle densities 
p(jt) = <5 77i <5(r — Tj) where the sum runs over all particles 

p L ( 1 r) = (p( 1 r)} L . (46) 



- n 

r 6 | : 



+ F SR (a a ,p b ,r a - r b ) 
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We keep the index L to remember that averages are taken for the finite- volume slabs 
(41). Hence expressing the sums in (45) as integrals on particle densities p(jt), the 
average force reads 

<F) L = -/dr fdv' SZ* Cl {t,t>) (47) 
Ja a (l) Ja b (l) |r - r I 

where c^(r, r') is the two-point charge correlation function 

c L (r, r') = (c(r)c(r')) L , c(r) = ]T e 7 p( 7 r). (48) 

7 

We now consider the average force by unit surface between two infinitely extended 
slabs at distance d by letting their transverse dimension L tend to infinity. We 
assume that the plasma phases are in fluid states homogeneous and isotropic in the 
y directions, namely the charge correlation has an infinite- volume limit of the form 

lim c L (r, r') = (c(r)c(r')) = c(x, x', |y - y'|). (49) 

L-+co 

For symmetry reasons, (F) L has no component parallel to the slabs and is directed 
along the perpendicular x axis. We therefore consider the x-component of the force 
per unit surface 

(/) ^ ft, ^ = ~h h {h V^f CL{X) y ' X ' 

r° f d+b f x-x' 

= - dx / dx> dy- c(x,x',\y\). (50) 

J -a Jd J [{x - x') 2 + |y| 2 J 

The last line results from the y translational invariance of the integrand in the limit 
L — > oo. Formula (50) remains valid if one replaces c(x,x', |y|) by the truncated 
charge-charge correlation function 

S(x,x',y) = (c(r)c(r'))-(c(T))(c(r')), r = (x,y), r' = (x',0) (51) 

with c(r) the microscopic charge density as in (7). Indeed, the y-Fourier transform 
of the Coulomb force reads 

/ dy C ' tk ' y [{x -^+\ yl 2 ? /2 = 27r ^ ~ x '>~ klX ' X ' 1 ( 52 ) 

and reduces to — 2tt when k = and x < x'. This implies that the charge density 
profile (c(r)) = c(x) does not contribute to the force because of the global neutrality 
of both plasmas 

fd+b 

dxc(x) = / dxc(x) = 0. (53) 

a J d 
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To take full advantage of the translational invariance in the y direction we rep- 
resent the y-integral in (50) in Fourier space : 

(/) = ^ Jdx Jdx'Jdke-^ x - x '^S(x,x',k), (54) 

where k = |k| and S{x,x',\s) = Jdke~ tk y S(x, x', y). The dependence of (/) = 
(/) (d) on the separation d between the two slabs occurs explicitly in (54) as well 
as implicitly in the charge correlation function S(x,x',k). The d dependence of 
the correlations between the two slabs A and B originates itself from the Coulomb 
interaction term Uab occurring in the total Gibbs thermal weight. 



3.2 The charge correlation between macroscopic conductors 

3.2.1 Mayer graphs and resummation of Coulomb divergences 

To determine the asymptotic behaviour of (/) (d) as d — > oo, we need to know that 
of the charge correlation between the sabs. For this we can use the technique of 
summation of Mayer graphs for Coulomb systems [29]. We recall that the two-point 
Ursell function, related to the densities p(i),p(J) and the two-particle density p(i,j) 

can be expanded in a formal power series of the densities by means of Mayer graphs. 
The basic Mayer bonds are 

f(ij)=e~^-l (56) 

where V(i,j) is the potential (43) and the weights at vertices are the densities 
p(i). Here i is a shorthand notation for the point (7^) in configuration space, and 
integration on configurations J^ 7 . f dri includes the summation on particle species. 
Diagrams have two root points % and j and m internal points which have to be 
integrated over. Each pair of points is linked by at most one /-bond and there are 
no articulation points 3 . Because of the long-range of the Coulomb interaction, the 
integrals occurring in every diagram diverge in the thermodynamic limit. The well 
known procedure to cure these divergences is to resum the chains built with pure 
Coulombic interaction bonds — /3e 7 ^e 7j .f (r^ — r 3 -). The chain summation leads to the 
following integral equation which defines the screened (or Debye-Huckel) potential 

$(r, r') = v(r - r') — dri K 2 (ri) v(r - ri) $(ri, r') = $(r', r) (57) 



3 An articulation point, when removed, splits the diagram into two pieces, at least one of which 
is disconnected from the root points. 
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which in turn is equivalent to the differential equation 



A$(r,r') -K 2 (r)$(r,r') 



47r5(r - r') . 



(58) 



In (57) and (58) 



1/2 



K 



(r):= 47r/?^e>( 7 r) 



(59) 
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can be interpreted as the local inverse Debye screening length ^(r) in the inho- 
mogeneous system. Once the screened potential has been determined, the original 
Mayer series for the Ursell function is reorganized into a series of prototype graphs 
which involve the two integrable bonds 



F = exp[-/3e 7i e 7j .$(r;, rj ) - Pv S r(^, 7j, \*i ~ ~ 1 + Pe^e^(r h Vj). 



The rules for prototype graphs are the same as for the original Mayer graphs except 
for the excluded convolution rule, namely the convolution of two F(i,j) bounds 
is forbidden (to avoid double counting of original Mayer graphs). Vertices receive 
density weight p( 7 r). In fact the particle densities are not known at this point, but 
they can be found self-consistently in principle from the first equation of the BGY 
hierarchy which links the densities to the two point correlation functions. Finally, 
the charge-charge correlation function (9) is related to the Ursell function by 



S(v, r') = <VyA*7 r)p(Y r>( 7 r, i r') + 5(v - r') ^ e>( 7 r) . (62) 



The second term in the r.h.s. of (62) is the contribution of coincident points. Tak- 
ing into account the sole bond F(i,j) (60) defines the Debye-Huckel or mean field 
approximation. 

In our model the particle densities p( 7 , x) as well as k(x) do not depend on the 
variable y. Written in Fourier space with respect to the y directions, the integral 
equation (57) and equivalent differential equation (58) become 



F{i,j) = -/3e 7i e 7i $(ri, ri 



(60) 



(61) 
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x f , k) = v(x, x', k) — — J dx 1 K 2 (xi)v(x,x 
' d 2 1 

— ^-k 2 -K 2 (x) x',k) — -4tt5(x - x'), 




dx 



(64) 



with 




(65) 
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Since the particles densities vanish outside the slabs one has 

k(x) = 0, if x < —a, < x < d, or x > d + b . (66) 

At the boundaries x = —a, 0, d, d + b, $(x, x', k) and its x- derivative must be 
continuous and <£>(x, x',k) — > 0, x — > ±oo. 

Solving (63) by iteration yields a series that can be shown to be convergent in the 
weak coupling regime. In particular <£>(x, x', k) has a bound uniform with respect to 
k and d, so that 

lim$(x, x',k) = $(x, x',0) < oo . (67) 

This means that the screened potential is short ranged (integrable) in the y direc- 
tion 4 . 

When d = oo, the system reduces to the single plasma slab A^. Its screened 
potential 

$° (x,x',k) = lim $(x,x',k), x,x' G A A (68) 

d^oo 

satisfies of course the integral equation 

$° A (x, x', k) = v(x, x', k) - / dxi«) 2 (xi)v(x, xi, k)$° (xi, x', k) (69) 

where k^(x) = lim^ 00 /t(x), x e A^, is the inverse screening length in the single 
plasma slab without any other electrical influence. 

3.2.2 Electroneutrality sum rules 

An important relation that holds in great generality in conducting phases is the 
electroneutrality sum rule [30] . It states that the total charge of the screening cloud 
around a specified charge ey at r' compensates it exactly, in terms of the Ursell 
function, it reads 

dr e 7 p(7, r)/i(7 r, 7 r) = -ey (70) 

7 

where the integrand in the left hand side is the charge density at r conditionned 
by the presence of a charge ey at r'. From (62) one immediatly sees that the 
charge-charge correlation obeys the sum rule 



J drS(r,r') = J dx J dyS(x,x',y) = (71) 



4 For finite d, $(x,x',y) decays faster than any power of |y|. For d = 00 the system reduces 
to a single plasma A with a hard wall without electric properties, for which $(x,x',y) decays as 
|y| -3 in the y-plane 
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or, in Fourier space with respect to y 

J dxS(x, x', k = 0) = . (72) 

An elementary derivation of (72) in the mean field approximation can be obtained 
from the differential equation (64). Integrating (64) on x leads to 

/•oo 2/ \ u2 r 

dx -^Hx, x',k) = l-^- Jdx $(x, x', k) (73) 

and in particular, for k = 

dx^-pt<5>(x,x',k = 0) = 1 . (74) 

J-oc 47r 

In the Debye-Hiickel approximation the Ursell function is given by the sole link (60) 
/i DH (7,x, 7 / ,a; / ,k) = -/3e 7 ey$(x, x', k) and from (62) 

S DK (x, x', k) = -P ^P- K ^ X 'h (x, x', k) + S(x - x')^P- (75) 
An An Air 

so that (74) implies (72) in this approximation. 

We will see that the electroneutrality sum rules will be at the microscopic origin 
of the universality of the Casimir effect. 

3.2.3 Asymptotic chain summation for large separation 

It is not easy to solve eq. (64) for two reasons: first the particle density profiles 
p(7, x) in the slabs are not explicitly known, and in any case the matching of solutions 
by continuity in the boundaries of different regions leads to cumbersome algebra. 
We shall instead follow another route by performing the chain summation only for 
the asymptotic part of the screened potential as d — > oo. For this it is convenient to 
split the Coulomb potential V into parts according to the location of its arguments 



v(x, x', k) 



' vaa(x, k) , a < x, x' < 0, 

vab(x, x' — d, k) , a < x < < d < x' < d + b, 

vba(% — d, x', k) , a < x' < < d < x < d + b, 

^vbb(x — d, x' — d, k) , d < x,x' < d + b, 



(76) 



and the same decomposition for the screened potential $ 5 . We can think of 

— (3vaa, —fivAB = —/3vba, —/3vbb as an expanded set of Coulomb bonds that will 



5 The functions vab, "bb, $>bb ■ ■ ■ (depending on d) refer to the system of the two 

plasmas under mutual influence with the allocation of particles in plasma B measured by their 
distance from the boundary at d (i.e. from to b). 
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enter into the chain resummation 6 . According to (65) the v AB potential can be 
written as 



v AB (x,x',k) = 
ke -kd 



2ire 



-k\x—x'—d\ 



ke 



-kd 



k 



2tt 



2vr 



-k\x\ 



27T 



-k\x'\ 



2ir 



vaa(x, 0)v bb (0, x'), v aa (x, 0) = v AA (0, x) = v BB (0, x) . (77) 



We call —fiv AB a traversing bond and the chains constituting — (3$ab traversing 
chains. Traversing chains have necessarily an odd number of traversing bonds. Let 
$S be the sum of traversing chains having n traversing bonds. Then introducing 
(77) and summing all chains of AA (BB) bonds attached to x (x') one obtains 



< §) = ^aa(x, 0, 2)<& Bfl (0, x', *) . 



qe 



(78) 



Note that because of the factor v AA (x, 0) the A A chains have necessarily an extremity 
located at the boundary x = of plasma A (likewise for plasma B). Thus, the 
dominant part of <&^ B (x, x', |) behaves as (see (68)) 



^ B (x,x'^) ~ ^-$°(a;,0,0)$ B (0,x',0), oo . 



qe 



2%d 



(79) 



One can obtain $ AB (x, x', V) by first attaching a traversing Coulomb bond to each 

extremities of $ AB (x, x', |) and summing the AA (BB) chains in plasma A (B) as 
before 



2vrd 



x 



$aa(x,0,1) 

J dx ± J dx 2 v AA (0,x 1 )'^^^(x 1 ,x 2 ^)'^^v BB (x 2 ,0)^ $ BB (0,x',Z) 



ge" 



2ttc? 



47T 



X 



«l( g 2) 
47T 



^Bs(0,x 2 )$A J 4(a;2,0, §) 



*bb(0,x',§) 



(80) 



where the second equality comes from (78) and rearranging the factors. By (69) the 
first square bracket in (80) is equal to 



^(0,0 2)-$ M (0,0 2) 



2nd 



+ 0(1) 



(81) 



6 Cases when ir, a;' are outside of the plasmas do not need to be considered since the potential 
will always be multiplied there by vanishing density factors. 
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since the screened potential $aa(0, 0, k) is finite at k = (see (67)), and the same 
estimate holds for the second bracket 7 . This yields a factor 

\2-wd J \ q J 2nd 

and leads to 

$ fi(^'^S^X(M 1 0)$° B (0 ) i',0) ) rf^oo. (82) 

Continuing by induction one sees that ^^ +1 \x,x' , ^) receives a factor 3 §^je~ 2nq 
and summing all the traversing chains gives the final result 

3>ab(x, x', §) ~ A7Td 9 sinhq ^A^ °> °) $ °b( ' ^ °)> - °o ■ (83) 

The screened potential of the joint system factorizes, up to a factor, in the screened 
potentials of the individual plasmas with a charge located at their inner boundary 

3.2.4 The asymptotic force 

We can now insert this behaviour into the basic formula for the force (54) 

(/> {d) = ^ Jh Jh'fa exp ( Hl "/-"' ) s ( *,j - i, i) 

x (0 y, ))=-m, ^oc. (84) 

To obtain the first equality we have changed the integration variables x' — > rr' — 
d, k — q/dm (54) and the second line follows from (75). In the last line we observe 
that both parentheses are equal to 1 because of the electroneutrality sum rule (74), 
whereas the q integral yields the value £(3) of the Riemann ^-function. 

This result deserves following comments. 

• The asymptotic value of the force has been calculated here within the Debye- 
Hiickel approximation using only the bond (60). It can be shown that the result 
(84) holds in full generality. In fact all the other Mayer graphs involving the 
bond (61) do not contribute to the asymptotics of the force. The factorization 
property (83) holds for the complete Ursell function as d — > oo in the form 

7 Strictly speaking, here k 2 (x),x e A a still depends on the influence of the second plasma As, 
but $,4,4(0, 0, k) tends to $^(0, 0, k) as d — > 00 and this does not modify the argument. 
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G A ( 1 ,x, la ,0,k = 0)G° B 



0) 



47r/3<isinhg 



e. 



e. 



(85) 



where G° A , G° B are correlations functions in the individual slabs A^, A B in- 
volving charges e 7a , e 76 at their inner boudaries that can be shown to satisfy 
the sum rule 



as a consequence of electroneutrality. 

• The result (84) exhibits universality in the sense that the asymptotic force 
does not depend on the microscopic composition of the plasma (chemical 
species, masses, charges) nor on the typical microscopic lengths (screening 
lenths £d, interparticle distances p -1 / 3 ). Moreover it does not depend either 
on the thickness a and b of the slabs: the Casimir force is entirely due to the 
charge fluctuations at their inner boundary. In this respect, if one neglects the 
microscopic lengths compared to their separation d, slabs of arbitrary thick- 
ness behave as infinitely thin conducting foils. Universality is a consequence 
of electroneutrality sum rules in the conductors. 

• If one compares (84) with (39) one sees that the extrapolation of Casimir 
calculation to the classical regime is larger by a factor 2 than what is obtained 
in the present classical microscopic model. The two approaches are based on 
different premises: (39) was derived from the electromagnetic field fluctuations 
but treating the metal as a macroscopic body without internal structure. The 
For (84) the force originates exclusively from the atomic fluctuations inside 
the metals, but since the dynamical degrees of freedom of the electromagnetic 
field have not been introduced, the force is purely electrostatic (longitudinal 
field) and the effects of the transverse components of the field (in particular 
the Lorentz force between fluctuating currents) are missing. This calls for a 
more complete model where quantum mechanics and photons are taken into 
account. 

3.3 Quantum corrections to the classical Casimir effect 
3.3.1 The complete model 

We consider the same two plasma slabs system A^, Ag as before, but now containing 
non relativistic quantum charges (electrons, ions, nuclei) with appropriate statistics 




(86) 
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immersed in a thermalized quantum electromagnetic field. The field is itself enclosed 
into a large box K with sides of length R, R^> L, a, b. The Hamiltonian of the total 
finite volume system reads in Gaussian units 

H = J2 ( p '~^ A(r<) ) 2 + J2 + E ^ WaUS ^' r *) + ^o ad • (87) 

27?v ^— ' rj — r,- ^— ' 

i 11 i<j J i 

The sums run on all particles with position r« and species index 7$; l /walls (7j, r«) is 
a steep external potential that confines a particle in or A#. It can eventually be 
taken infinitely steep at the wall's position implying Dirichlet boundary conditions 
(i.e. vanishing of particle wave functions at the faces of the slabs). 

The electromagnetic field is written in the Coulomb (or transverse) gauge so that 
the vector potential A(r) is divergence free and ifg ad is the Hamiltonian of the free 
radiation field. For it we impose periodic boundary conditions on the faces of the 
big box K 8 . Hence expanding A(r) in the plane waves modes k = (^rS ^r 2 -, 
gives 

A(r) = (^) 1/2 ^ 9 (k)^(a ; , A e-'- + a k , A e"") (88) 

k, A 

H™ d = E^< A a k , A . (89) 

k,A 

In (88), <7(k), g(0) — 1, is a form factor (ultraviolet cut-off) needed to make sense of 
the Hamiltonian (87) 9 . In our case we are concerned with the asymptotics d — > oo 
which is related to the k — > behaviour, hence the final result is expected to be 
independent of this cut-off function. 

We suppose that the matter in the slabs is in thermal equilibrium with the 
radiation field and therefore introduce the finite volume free energy of the full system 
at temperature T 

$*,L,d = -£; B TmTre-^ (90) 

where the trace is carried over the particles' and the field's degrees of freedom. The 
force between the slabs by unit surface is now defined by 



with 



f(d) = lira lim f R , L (d) (91) 



fnAd) = • (92) 



8 Periodic conditions are convenient here. We could as well choose metallic boundary conditions 
as in chapter 1. Since the field region K will be extended over all space right away the choice of 
conditions on the boundaries of K are expected to make no differences for the particles confined 
in the slabs. 

9 Removing this ultraviolet cut-off is still an open problem. 
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Adding and substracting the free energy of the free photon field in (90) leads to 



®r,la = -k B Tln 



Tre 



-0H 



k B T In Z™ d 



(93) 



where Zg ad is the partition function of the free photon field in the volume K. Since 
the last term of (93) is independent of d, it does not contribute to the force (92). 
Therefore one has 



1 d 

f(d) = k B T lim lim y^^m 



Tie-? 11 

^rad 



(94) 



The thermodynamic limit is expected to exist when one at least of the species is 
fermionic (say the electrons) without regularization of the Coulomb potential at the 
origin. 

The situation is very similar to that considered by Casimir: The field extends 
over all sides of the two conducting plates, but the latter are now described at 
the microscopic level including all the particle- field interactions. At this point, the 
following observation is important. The Bohr-van Leeuven theorem [31] states that 
classical matter in thermal equilibrium decouples from the transverse part of the 
(classical) electromagetic field. We recall the central argument which is very simple. 
In the phase space integral of the classical Gibbs weight 



dr / dpexp 



-P 



(P ~ | A(r) 

2m 



exp(-/3C/(r)) 



(95) 



one can perform the momentum integration first, shifting the variable p to p — 
§A(r) fot fixed r so making the integral independent of the vector potential A. 
It is therefore expected that in the classical limit the thermal statistical averages 
calculated with the QED Hamiltionan (87) reduce to those obtained by the purely 
classical model of section 3.1 based on the sole Coulomb electrostatic interaction. It 
is thus of interest to study the Casimir force from the full quantum model (87) in 
the semi-classical (or high temperature) regime. 

In addition to the photon thermal wave length /3Hc, quantum mechanics intro- 



1/2 



duces particles' de Broglie thermal wave lengths A 7 = h . The semi-classical 

regime is the situation where all thermal wave length are much smaller than d and 
the slab thicknesses a, b and the parameter a (21) is also small. We establish be- 
low that, in this regime, the quantum corrections to the classical Casimir effect are 
small. More precisely, they do not contribute to the dominant order d~ 3 so that one 
can write 



fid) 
R((3, h, d) 



C(3) 

87^3 

e>(cr 4 ) . 



+ R(p, h, d) 



(96) 
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The remainder R({3, H, d) includes the quantum effects, but it cannot modify the 
amplitude of the d~ 3 term to match the high temperature result — 4 ^] 3 (39) obtained 
by Casimir's method ignoring the microscopic fluctuations inside the conductors. 

3.3.2 The path integral representation and the effective electric and 
magnetic potentials 

For the sake of simplicity, we shall sketch the procedure with the following specifi- 
cations (for more details, see [32]). 

• The spins of the particles are ignored (as it is already the case in the Hamito- 
nian (87)). They could be introduced but we expect that their contributions 
will not modify the asymptotic force. 



• The Fermi or Bose statistics of the particles will not be taken into account. 
The use of Boltzmann statistics is in order in the non degenerate regime that 
we are considering, but requires the presence of a short range repulsive po- 
tential t>sR(7i;7j) \ r i ~ r j\) to assure stability as in the classical case. Anyway, 
exchange effects across the two slabs will be negligible at large separation. 



• We treat the electromagnetic field classically, which is justified by (3hc <C d 
(i.e. a < 1, see (21)). This amounts to replace the photon creation and 
annihilation operators by complex numbers a£ A , «k,A- 

The formalism adapted to the investigation of the high temperature (or semi- 
classical regime) is the Feynman-Kac-Ito path integral representation of the Gibbs 
weight. For a single particle in an external potential V^ ext (r) the formula for a 
configurational diagonal matrix element is 




3/2 



exp —f3 



2nX 2 J 
i 



J dsV™\{v + \i{s)) - A(r + A£(*))-de( 



s) 



(97) 



Here £(s), < s < 1, £(0) = £(1) = 0, is a closed dimensionless Brownian path 
and D(£) is the corresponding conditionnal Wiener measure normalized to 1. It is 

Gaussian, formally written as exp ^— \ ^ ris^( s )> w ^ n zero me an and 



covariance 



/ 



(98) 
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where ^(s) are the Cartesian coordinates of £(s). 

In this representation a quantum point charge looks like a classical-like structured 
charge at r with an internal degree of freedom, the random charged filament £(s) 
whose extension is given by the de Broglie length A (the quantum fluctuation). The 
magnetic phase in the bracket of (97) is a stochastic line integral: it is the flux of 
the magnetic field across the closed filament. 

This is readily generalized to a system of n interacting particles: the two first 
terms of (87) yield the following expression of the Gibbs weight in the space of 
filaments 



exp 



i<j j=l 3 J ° 



where 



Vc^iAi^jAj) = [ ds- 



(99) 
(100) 



|ri + A 7i ^(s) -Tj - A 7 .^(s)| 

is the Coulomb potential between two filaments. 

A remarkable fact about the representation (99) is that the exponent is linear 
in A and its Fourier amplitudes (contrary to the Hamiltonian (87) written in op- 
eratorial form). Since the statistical weight e~@ H ° a is a Gaussian function of these 
Fourier amplitudes, it makes it possible to perform the partial trace over the field 
degrees of freedom in (94) by the fact that the Fourier transform of a Gaussian is 
a Gaussian, namely for a n x n positive definite hermitian matrix C and complex 
vectors z = {zi}, J = {Ji}, i — 1, . . . , n 

J l\ n Det(C) e • tiUij 

This formula is applied to the calculation of this partial trace in the form of the 
normalized Gaussian average 

< - i^infv^ - ■ (io2) 



The result is 



f[ e-^'f™™^' ^ \ e -^<l e ^j W ™( r i£i> r i£j) . (103) 
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In (103) W m is a double stochastic integral 

1 



3> 



m li m lj c 



dk 

(2^ 



(104) 



where 



GHk) = ^t^C(k), C(k) = ^-^ (105) 



is the free field covariance and 5^(k) the transverse Kroneker symbol. Summation 
on the Cartesian componant indices /i, v = 1,2,3 is understood. We see from (103) 
that W m can be interpreted as a pairwise effective magnetic interaction between the 
random filaments £ { and ^ mediated by the vector potential A 10 . This interpreta- 
tion becomes manifest when we look at the long distance behaviour of W m , 

W m (r 1 ,£ 1 ,r 2 ,€ 2 ) ~ — : 1 - I ^i(si) • / d£ 2 (s 2 ) 

(AiCi(si) • V ri ) (A 2 £ 2 ( S2 ) • V r2 ) ----- 



ri - r 2 



(106) 

as |ri — r 2 | — > oo. This is obtained by keeping the most singular term in the 
small k expansion of the exponentials in the integrand in (104). (Here we have 
omitted additional 0(r~ 3 ) terms coming from the transversality condition in (105).) 
Since filaments are closed one has d£(s) = 0, and this implies that the first non 
vanishing term is bilinear in £ : and £ 2 . If one interprets formally 



f^(x-r-A^))v( S ), v( 
Jo 



j(x) = e / ds6(x - r - A£( S ))v(s), v(a) = (107) 



as the current density carried by a wire at r, one sees that (106) has (up to a factor) 
precisely the form of the classical magnetic energy of such a pair of current wires. 
A comparison of (100) and (104) shows that the ratio W m /V c ~ (Pmc 2 )' 1 is a small 
quantity at high temperature so that purely electrostatic effects will be dominating. 

Having now identified the basic effective pair interactions between the filaments, 
namely the electrostatic potential (100), the magnetic potential (104) and a possible 
non electric short range potential as in (43), it is possible to proceed as in the classical 
treatment of section 3.2.1. Indeed if one considers the auxiliary system of filaments 
r, £ described in an enlarged phase space of classical-like particles equipped with 



10 



The product in (103) contains the magnetic self energies of the filaments. 
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an internal degree of freedom £ and the above pair interactions, all concepts and 
methods of classical statistical mechanics apply. In particular we have the density 
p(7, r, £) of filaments of species 7 and the filament-Ursell function /i(7, r, 7', r', £') 
which can be expanded in Mayer graphs. The Coulomb part (100) still decays as 
r _1 so that graphs still suffer from the long range divergences, and chain summation 
have first to be performed. There is an important observation to be made at this 
point: from the Feynman-Kac formula the potential (100) inherits the equal time 
constraint, i.e. every element of charge eiAid^Si) of the first filament does not 
interact with every element e 2 \ 2 d£ 2 (s 2 ) as would be the case in classical physics, 
but only if si — s 2 . It is therefore of interest to split he coulomb potential into 

V(i,j) = V aec (i,j) + W c (i,j) (108) 

where 

V c i cc (i,j) = ds! ds 2 - (109) 

Jo Jo \*i + K£i{ s i)- T 3-Kwv\ 

is a genuine classical electrostatic potential between two charged wires and 

W e (iJ)= [ 1 ds 1 f ds2 (5( Sl - s 2 ) - 1) - 1 (110) 

is the part of V(i,j) due to intrinsic quantum fluctuations (W c (i,j) vanishes if h is 
set equal to zero). Because of the identities 

f ds 1 (S(s 1 -s 2 ))-l)= f ds 2 (5( Sl -s 2 ))-l) = (111) 
Jo Jo 

its large distance behaviour originates again from the term bilinear in £ : and £ 2 in 
the multipolar expansion of the Coulomb potential in (110) 

~ [ d Sl [ ds 2 (5( Sl - s 2 )) - 1) (A^( Sl ) • V ri ) (\ 2 £(s 2 ) • V ra ) , 1 , ■ 

JO JO l r l — r 2| 

(112) 



It is dipolar and formally similar to that of two electrical dipoles of sizes ei\i£ 1 
and e 2 X 2 $ i 2- The chain summation of the Coulombic part V c \ ec (i,j) will lead to a 
screened potential $ e iec(^i)- Then, as in the classical case, one can introduce the 
Mayer diagrammatics in prototype graphs with bonds 

F{i,j) = -/3e 7i e 7 .$ e iec(i, j) (113) 
F R (i, j) = exp[-/5e 7i e 7i ($ c i cc (i, j) + W c (i,j) + W m (i,j))] - 1 + /3e 7i e 7 .$ e i ec (i, j) . 

(114) 
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Weights at vertices are the filament densities p(7, r, £) and one has the excluded 
convolution rule for the bonds F(i,j). This is similar to (60) and (61) with the 
additionnal occurence in F R (i,j) of the electric and magnetic multipolar interaction 
W c and W m due to intrinsic quantum fluctuations of the filaments. 

The rest of the analysis is as follows. The potential $ e i ec is a classical mean field 
potential for structured charges. Its large d contribution can be extracted exactly as 
in sections 3.2.3 and 3.2.4 leading to the same universal result — 8 ^] 3 as (84) (using 
that electroneutrality sum rules also hold in the system of filaments). 

In contrast to the classical situation, at large distance the bond F R has now a 
dipolar behaviour F R (i,j) ~ — f3e li e lj {W c (i, j) + W m (i, j)) ~ |rj — rj| -3 . Concerning 
the large separation d, this bond (in the partial Fourier representation with k = q/d) 
proves to decay as 0(l/d), exactly like $ e i ec (q/rf). This is easily seen for W c (q/d): 
using (65) with k = q/d and expanding for large d, terms of order 0(d) and 0(1) 
identically vanish because of (111). Basically the same holds for W m (q/d) except 
that the vanishing of these dominant terms follows from the vanishing of the line 
integral over a closed filament f Q d£(s) = 0. Therefore, the dominant behaviour of 
the complete filament Ursell function Hab{1i x, l\ x', ^) contains in addition 
to (85) (expressed with filament degrees of freedom) contributions built with sole 
traversing links W Cj ab(1j 2) and W mj AB(l,2). To point 1 in A one can attach all 
possible prototype graphs of system A and likewise to point 2 in system B resulting 
in the additional terms 



However, these new contributions have no effect at leading order 0(1/ d 3 ) when 
integrated in the force because of the electroneutrality sum rule in the form (71). 

The major difference with the term (85) containing the correlation functions 
G° A , G° B and leading to the asymptotic result (84) is that the latter functions are 
constrained by the excluded convolution rule as they are attached to the traversing 
bond Fab, and that they satisfy in turn the sum rule in the form (86), yielding 
universality (instead of the vanishing) of the 0(d~ 3 ) contribution. 

3.3.3 Concluding remarks 

We have shown that quantum corrections to the classical Casimir effect can con- 
veniently be worked out from the electric potential W c (i,j) (110) and magnetic 
potential W m (i,j) (104) occuring naturally in the functional integral representation 
of the system. These potentials account for the random charge and current inter- 
actions generated by the quantum fluctuations. The size of these interactions at 
long distance are measured in terms of the thermal wave length A of the particles. 
It is remarkable that the dominant d~~ 3 term is exactly the universal classical one, 




dl J d2 [p A (l)h° A (a, 1, k = 0) + 5(a, 1)] (-/3e 71 e 72 ) 
(W c + W m ) (1, 2, a) [p° B (2)h B (2, b, k = 0) + 6(2, b)} 
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still independent of the microscopic details of the conducting phase (interparticle 
distance, screening lengths, thermal de Boglie lengths). This is again a consequence 
of the electoneutrality sum rules extended to the quantum sytem. This calculation 
strongly confirms that (84) is the correct value of the Casimir force at high temper- 
ature. One must conclude that the discrepancy with (39) is not due to the omission 
of the transverse part of the electromagnetic interaction in the classical Coulombic 
models, but should be attributed to the very fact that fluctuations inside the conduc- 
tors are ignored in the calculation leading to (39) 11 . In other words, the description 
of conductors by mere macroscopic boundary conditions is physically inappropriate 
whenever the effect of thermal fluctuations on the force are considered. 

One the other hand, recent experiments validate the zero temperature formula 
(20). In [25] the authors find an experimental agreement with the value of Casimir 
force's strength ir 2 hc/240 to a 15% precision level. This indicates that fluctuations 
in conductors are drastically reduced as the temperature tends to zero and possibly 
have no more effect on the force at T = 0. In (96), quantum effects appear at the 
subdominant order d~ A . One may imagine the following scenario: as the temperature 
is reduced, the classical term (~ T/d~ 3 ) decreases whereas the term R((3, h, d) ~ gT 4 
approches the Casimir vacuum value — (20). Understanding the cross over from 
(96) to the zero temperature formula of Casimir is an open problem. 



4 Dispersion forces 

4.1 Van der Waals - London forces in vacuum 

It was a great achievement of the early quantum mechanics to show that there is a 
general force of attraction between atoms or molecules even if neither has a perma- 
nent average dipole moment. The force comes from the intrinsic quantum fluctua- 
tions of the charges inside the atoms, distributed according to the wave function of 
the atomic states. For this reason such microscopic forces can also be considered to 
be of Casimir type. They are called dispersion forces because they involve the po- 
larizability of the atoms which is also related to the refractive index of the medium. 
London's derivation is part of the standard education in quantum mechanics [34] 
and we just recall its principle. 

For instance, for two Hydrogen atoms in their ground state and ip^ with 
infinitely heavy nuclei located in and the mutual Coulomb interaction behaves 
as a dipolar potential at large atomic separation r=|r|^oo, r = — r B , 

p\at—at 

V at - at (r) ~ , D at ~ at = e 2 [y (A) • y (B) - 3(y (A) • r)(y (B) • r)] (115) 

where y^ and are the relative electronic coordinates and r = r/r . The van der 
Waals potential uwij") obtained by treating V at ~ at as a second-order perturbation 

11 In fact it is well known that there are long range field correlations inside the conducting media 
when T^O, see [33] for electrostatic potential fluctuations and [32] for transverse field fluctuations. 



31 



reads 

uw(r) = -C w /r fi 

° W = ^ E + E -2E n >0 - (116) 

In (116) E m , m = 0, 1, . . ., are the eigen energies of the Hydrogen atom (repeated 
according to their multiplicities) and the sum runs on all excited states ipffl <E> 
■ipn i (m,n) 7^ (0,0), of the two atoms (the notation includes the integral on the 
continuous part of the spectrum). 

It will be usefull in the sequel to relate the van der Waals amplitude Cw to the 
polarizability a(u) of individual atoms. The latter is defined as the linear response 
of the electric dipole ey of the Hydrogen atom in its ground state to an applied 
oscillating electric field E e luJt : 

(Mt)\ey\Mt)) ~ a(w)E e <wt , E - (117) 
and has the value [35], sec. 98, [36], sec. 92 

a{u) ~ 2e h (^-w-(m 2 • (118) 

Taking the 3-axis along r, one has D at ~ at = e 2 {y^ y^ +y^ y^ — 2y^y^), explic- 
itating the matrix elements in (116) and exploiting that ipo is spherically symmetric, 
one can write Cw in the form 12 

c =6 4 V- I(^0|^l^)| 2 |(^0b3 g |^)| 2 fn9 x 

W jka, E m + E n -2E ' 1 ) 

The connexion between (118) and (119) is provided by an application of the identity 

1 2ab r° 1 
V+b ~~J du (a 2 + u 2 )(b 2 + u 2 ) { ' 

yielding the desired relation 



3h f°° 

Cw = — / du atA(iu)aB(iu) . (121) 

71" Jo 



12 Because of rotation invariancc, (118) and (119) could be written in terms of any Cartesian 
component of y. 
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4.2 Van der Waals - London forces at finite temperature 

The above calculation disregards all many-body and temperature effects which are 
present when atoms are in a thermal equilibrium state at temperature T and density 
p. It is therefore of interest to study the effective atom-atom potential in a fluid 
with non vanishing T and p. Two questions arise: does this effective potential still 
decay as and if it is the case, what are the temperature and density corrections 
to van der Waals amplitude Cw (116) for an atom pair in empty space ? The first 
question is in fact not trivial as illustrated by the following simple and apparently 
sensible reasoning which turns out to be incorrect. At non zero temperature the gas 
always contains a fraction of ionized free charges that will be the source of screening. 
As commonly done, one could take the screening effects due to these free charges in 
the medium into account by replacing the bare Coulomb potential between charges 
of different atoms by the screened potential obtained in the usual Debye-Huckel or 
RPA mean field theory. Since the latter decays exponentially at large distance, the 
effective attractive interaction between two atoms would also be exponential. This 
reasonning predicts that the 1/r 6 van der Waals forces should disappear as soon 
as there is a fraction of thermally ionized charges, a false conclusion. What really 
happens for the quantum gas becomes particularly clear when we use the Feynman- 
Kac path representation of the correlations presented in section 3.3 (here we consider 
the pure electrostatic model without the photon field). In this formalism point 
particles appear as random charged filaments carrying mutipoles moments. Then 
the clue is given by the decomposition (108) of the Coulomb interaction between 
filaments into a classical part (109) and the proper quantum fluctuation part W c 
(110) not reducible to any kind of classical behaviour. As (112) shows W c decays as 
a dipolar interaction r~ 3 and is the source of an algebraic decay of the correlations 
between atoms. In fact the bulk particle correlations decay as r~ 6 (because of 
rotation invariance, W c does not contribute and the dominant part of the decay 
at long distance is determined by ~ r~ 6 ) (see [37] and [18] for a review and 
references therein). 

It is interesting to discuss the status of the van der Waals forces in the Saha 
regime when atoms and free ionized charges are in thermal equilibrium (equilibrium 
ionization phases). Such phases occur provided that the temperature is sufficiently 
low to prevent full ionization an the density is low enough to have non overlapping 
atomic wave functions. For instance in the partially recombined Hydrogen plasma 
the densities of ionized electrons (e) protons (p) and the density of Hydrogen atoms 
in their ground states (at) are characterized by the expression that they would have 
for ideal gases 

p ' = p ' = j2^x^r^" = < } t < 122 > 

with A e , A p the thermal wave lengths of the electron and of the proton. These 
densities have to be equal because of neutrality and are denoted p l f , the density of 
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free charges. The ideal atomic density is 



ps = wkw^ lE °~ 2 *' = ^ 3 ' 2 [""] 2 e ' 0E " (123) 

with Eq < the ground state energy of the Hydrogen atom and A the thermal 
wavelength associated with its reduced mass. The chemical potential \x determines 
the total average particle number. If one sets 

fi = ntf) =E + k B T In w (124) 

with w a fixed parameter < w < oo, one sees on (122)-(123) that the ideal 
densities of protons and electrons become of the same order as the atomic density. 
The system behaves as a mixture of protons, electrons and Hydrogen atoms in their 
ground state. This describes the Saha regime of equilibrium ionization, with w fixing 
the relative proportion of free charges and atoms, 



rfS=«('=^)'V (125) 
V. y/m e m p J 

We report here without proofs the main result of the work [21]. We consider 
the proton-proton correlation p pp (r) at fixed temperature and total density p (r is 
the distance between the two protons). As said before in the fluid phase and for all 
positive values of T and p, Ppp T (r) behaves asymptotically at large distance as [37] 



We show that in the ionization equilibrium regime the coefficient C (T, p) takes the 
form of a sum of three contributions 

C(T,p) = {[pf] 2 Cf-f(T)+pfp%Cf- at (T) + [f%] 2 C«-*(T)} (1 + 0(e-°"« T )) 

(127) 

giving rise to three large-distance effective potentials 

(~if—f (if— at (iat— at 

uf-f(r) = - — , „/-( r ) = -_ u — ( r ) = __. (128) 

The three terms reflect the fact that, in the limit, a proton can be thought of as 
either being free or belonging to a Hydrogen atom. At lowest order in density the co- 
efficient is quadratic in the ideal free and atomic densities up to exponentially small 
terms (9(exp(— c/ksT)), c > 0, that include all higher-density effects. The factors 
C*~*(T), C^~ at (T) and c at ~ at (T) represent the effective interaction strengths be- 
tween two free protons, a free proton and an atom, and two atoms. They are still 
temperature dependent and have the asymptotic values as T — > 
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C f ~ at (T) 



h 2 



12 knT \ m 



p e/ m^O u=l 



x 



6k B T l2(k B T) 2 ,„ f . eM ' 



E m - E (E m - E ) 2 (E m - E ) 



(130) 



C at-at {T) = J- \^ a) ®4 b) \D at - at \^®^)\ 2 

1 2k B T 



X 



E m + E n — 2_E (-E m — E )(E n — E ) 



+ 0(exp(-5/fc B T)) 



(131) 



The results are the same for the electron-proton or electron-electron correlations: in 
the atomic limit pfj 1 (r) and p£} T {r) behave as Ppp T (r), with the same amplitude 
C(T,p). The temperature power-law corrections in (130) and (131) come from 
partial screening due to the presence of free electrons and protons. The exponentially 
decaying terms C(exp(— 5/k B T)) include the contributions of the excited states with 
< 5 < E x - E . 

The first noticeable point is that there is a van der Waals type effective potential 
uf~f(r) between unbound charges: quantum screening reduces the bare Coulomb 
potential r _1 to r~ 6 but not further because quantum charges together with their 
screening clouds behave as fluctuating dipoles. uf~f (r) has the remarkable property 
to be independent of the charge species (proton or electron). Thus it is attractive 
irrespective of the charge signs. 

The potential u f ~ at (r) with coefficient (130) results from the coupling between 
free-charge and atomic-dipole fluctuations. The quantity Djj~ at refers to the inter- 
action between the atomic dipole ey at the mass-center position with a reference 
dipole eu M located at r a . It is defined as in (115) by 

D f - at = e 2 [u(«) ■ yW - 3(fiW • r)(y<« • f)] (132) 

where for \i = 1, 2, 3 u M is a unit vector along the //-axis. In vacuum the charge 
dipole potential decays as r~ 4 . More precisely if a charge e is placed at large distance 
from an Hydrogen atom in its ground state in vacuum (Stark effect due to a localized 
charge) one finds 



V'-~(r)~-% Co-E '^f''? "' >0 - < 133 > 

m^O 



In the medium the decay changes from r" 4 to r~ 6 , becoming again of dipolar type. 
Finally the traditional van der Waals contribution (131) reduces to its vacuum value 
C\y as T — > 0. 
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Concerning the temperature dependence of the various amplitudes, a comparison 
with the van der Waals coefficient Cw = 0(1) gives 

C f - f (T) ~ {j§j C w, C^\T) ~ ^C w . (134) 

The larger size of the coefficients involving free charges has to be traced back to 
the larger size of their dipole fluctuations (of the order of A 2 ~ T _1 compared to 
atomic ones proportional to the square Bohr radius ~ a 2 B ). As a consequence, in 
the ionization equilibrium phase where p 1 ^ and p^t are °f the same order, the large- 
distance behaviour of the particle-particle correlation is dominated by the effective 
interaction between free charges. In the purely atomic phase (p l j exponentially 
smaller than p 1 ^), standard van der Waals forces are the only relevant interactions. 

We emphasize that these results are exact: they are derived in [21] by means of a 
renormalized virial expansion for the quantum Coulomb gas taking into account all 
the different effects stemming from the Coulomb potential at various scales (atomic 
binding, collective screening, polarization) in a systematic and coherent way. 



4.3 Retarded van der Waals forces in vacuum 

In the London calculation of van der Waals forces only the electrostatic part of the 
interaction (the Coulomb potential) has been taken into account. In fact one has to 
put at work the full electromagnetic interaction including the retardation effect due 
to the finiteness of the velocity of light. For this one has to add the transverse part 
of the field by writing the coupling of a charge to the vector potential as usual 



— (p-^A(r)) 2 = ^ + 



2 

(135) 



-£ P ' AW+ 2^ |A(r » 



Casimir and Polder [1] treated the second term of the r.h.s. by the standard methods 
of perturbation theory, and were obliged to go up to the fourth order since the 
result is expected to be proportional to the fourth power e 4 of the electronic charge. 
In terms of elementary processes it corresponds to Feynman diagrams where two 
photons are exchanged between the two atomic electrons. In order to interact with 
the field and return to its ground state each electron must necessarily emit or absorb 
two photons. After an elaborate calculation, Casimir and Polder found 

TA . . 23hc 

V(r) ~ -a A a B ^ , (136) 

where a:^, are the static polarizabilities of the atoms (each of them proportional 
to e 2 ). 

Following [4], sec. 3.11, we shall present another derivation, inspired from 
Casimir's second paper on the subject [3] which gives more physical insight. Here 
van der Waals forces at zero temperature are also seen as arising from the vacuum 
fluctuations of the field. The model is as follows. 
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• One treats the atom or colloidal particle at as globally characterized by 
a classical electric dipole u^it) = i>A^A{t) of orientation v a, \v a\ = 1, size 
VA(t), and by a frequency dependent electric susceptibility a^(a;). The second 
atom at r B is described likewise (the microscopic electronic structure of these 
particles are not spelled out). 

• These dipoles are fluctuating. They are induced as linear response of the atoms 
to the vacuum fluctuations of a free quantum electromagnetic field. 

• Each dipole acts as a source of a radiation field. The total field is considered 
to be the sum of a free quantum field (vacuum field ) and the field due to the 
sources. 

• The energy of the two atoms at distance r = \ta — i\b| is calculated as the 
energy of the second dipole i>b at in the field due to the first dipole v a at 
Va plus the vacuum field. 

The free electric field with periodic boundary conditions is obtained from (88) 
EoM) = -^AoM) = E[ E Jk,AW^ kt + E o, k ,A(r)e^ kt ] 

k,A 

Ej M (r, t) = -i (^< A (r, t)e k (A)e^ = (E^r, t)f (137) 

where E + (E~) denotes the positive (negative) frequency amplitudes of the field. 
The induced dipole of atom at by the k mode of this field at linear order is 

= < k ,A^ kt + ^, M e-^ 
< k ,A = *EJ M W = W (138) 

so that v\ k A e ltJk * represents an oscillating dipole of frequency cj k , orientation i> A = 
e k (A) and strength z/+ kA = a A (uj k )E+ MX (r A ). 

At this point we can recall the classical formula for the positive frequency am- 
plitude of the electric field at r radiated by an oscillating dipole v(t) located at the 
origin, [38], sec. 9.2, 

E(r) = E far (r) +E ncar (r), r = |r| 
E k J(r) = [3(u -f)f - ,]e + e -*r (139) 
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This formula is valid when the typical wave length A = c/2huj of the radiated wave 
is much larger than the atom size 13 . For r — > oo, E far correspond to a spherical 
outgoing wave whereas for r small E near reduces to the electrostatic dipolar field. 

Calling E j4 (r) the field (139) due to the dipole (138) located at r A , the total field 
at r is then 

E(r) = E (r) + E A (r) . (140) 
This field in turn induces a dipole 

< M = »bMK,x(*b) (141) 

on atom B at r^. The energy of the dipole pair is obtained in bringing the second 
atom B from infinity to in the field (140), keeping the first atom fixed at r^. 
This energy is — \v B • E(r#), r = \r A — r B \, see [38], sec. 4.7, 14 so that with (141) 
the contribution to the energy of the k, A mode is — ^aB(<^k)|Ek,A(i'B)| 2 - Finally, to 
obtain the total average interaction energy V(r), we have to sum over all modes and 
take the mean value in the vacuum state |0) of the photon field 

V (r) = -\j2 a BM(0\ |E k , A M| 2 |0> . (142) 

k,A 

Since the vacuum |0) is invariant under the free field time evolution, V(r) is time 
independent. In fact, we have already set t — along its derivation. 

Keeping only the terms that are bilinear in the polarizabilities gives 15 , 

V(r) = -\Y. aB ^) (°l [ E oM • E A (r B ) + E A (r B ) • E (r B )] |0) (143) 

k,A 

From now on to find the r dependence of V(r) is a matter of calculation (see [4]). 
The explicit formulae for the fields are introduced from (137), (138) and (139). Since 
the expression to be averaged is quadratic in the photon creation and annihilation 
operators, the vacuum expectation is easily computed from (0\ak t \a^ x \0) = 1. One 
finds in the infinite volume limit and performing the polarization sums 

V(r) = -j^Rej J™ dk k'>uj k a A (u> k )aB(uJ k )e-< k '' j tia ie ik -" 

[w^s + 1 1 - 3 * ' f > 2 i iwr + W?) } ■ (144) 

13 In our situation, the dipole is an operator proportional to a£.A- Since the Maxwell equations 
governing the quantum field (in the Heisenberg picture) depend linearly on the amplitude of the 
source, the solution for the radiated field is the same as in the classical case. 

14 The factor 1/2 is due to the fact that the dipole vb is induced and not permanent. 

15 To keep the term lE^r^)! 2 quadratic in a^Wk) would not be consistent with the linear 
response assumption. 
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The angular integration can be performed with the final result ( a>k = ck) 
V(r) = ~Rey^duu 6 a A (u)a B (u)D(<f)Y 
= e - (__ + _ + _____) 

In order to find the asymptotic behaviour of V(r) as r — > oo one needs to specify the 
u; dependence of the polarizability. To keep the discussion simple consider a single 
term of the expression (118) 

, x O ^ 4 (^n-^o)|(^o|l/3|^)| 2 ,^ 

a(uj) ~ 2e — — —- — — rr— (14b) 

assuming that a particular atomic transition between the two states 0, n gives a 
dominant contribution to the polarizability. One notes that a(u) has no poles in 
the complex uj plane and its value is real on the imaginary axis uj = iu. Then by 
a 7r/2 rotation of the integration line (0, oo) we can carry the integration in (145) 
along the imaginary axis uj — iu giving 16 

V(r) = -^Re ji duu e a A (iu)a B (iu)D(if)^ . (147) 

• Long distance 

Since the function D {^f-) has a factor exp(— ^f) it is clear that only the static 
susceptiblities a = a(0) will contribute as soon asr> hc/\E m — E \ . Indeed 
changing the variable ^ = v one has 

V(r) ~ -^-o^a B Re \i [ dvv 6 D(iv) 1 — — ^ ^ tt^ag, r — > oo (148) 
7rr' I Jo J 47T r ' 

which is the result of the Casimir and Polder calculation. 

• Short distance 

For small r the dominant contribution comes from the x~ 6 term in (145). This 
gives in (147) 

3h f°° _ 2w 

V(r) ~ g / dua>A(iu)a>B(iu)e c 

7rr Jo 

/•oo 

/ duaA(iu)aB(iu) . (149) 

•/ 



3/1 



7rr 6 



The last line holds when r <C hc/\E m — E$\. This is London's result (121): it 
can be formally obtained from (147) by taking the static limit c — > oo where 
only the contribution of the near field term in (139) survives. 



16 



The large quarter of circle closing the positive quadrant does not contribute. 
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4.4 Forces between dielectric bodies 

This topic is the subject of a very large literature and we will only summarize some 
main lines of the theory. 

4.4.1 Dilute dielectric bodies at zero temperature 

One can easily obtain an information on the type of decay between two dielectric 
slabs containing polarizable atoms at very low density. If correlations between atoms 
are ignored one can simply sum all the van der Waals potentials between atoms pairs. 
Let pa, Pb be the atom densities in the two slabs and fore sake of generality, assume 
that the pair potential behaves as — at large distance for some rj > 4, B > 0. 
Then the total potential energy per unit surface is 

u(d) ~ - lim f dv [ dr> 1 



l^oo L 2 J Aa J Ap |r-r'h 



V A </ 1\ B 

r-0 /-CO 



= -Bp A p B I dx i dx [ dy— 

J-oo Jd J {{X-X') 2 

= -B PAPB — JjxJ d dx^-^ 
2-kBpaPb 1 



+ y 2 y>/ 



2 



(77 — 2)(t7 — 3)(?7 — 4) rf?- 4 



(150) 



The force is f(d) = — and we find from (116) the short distance (non retarded) 
behaviour (77 = 6) 

tfA\ kCwPaPb 1 , . 

M — T s (151) 

and from (148) the long distance (retarded) behaviour (rj = 7) 

f fj\ 23hcp A aApBaB n - , 
fid) m , . (152) 

Whereas the powers of decay laws are correct, it was experimentally recognized in 
the early 1950s that the amplitudes found from (151) and were (152) not. This is not 
astonishing since this calculation is performed at lowest order in the atomic densities 
disregarding all correlations between atoms. This motivated Lifshitz in 1956, [39] 
and [40] sec. 90 (and many authors after him) to develop a theory taking into 
account the physical properties of the dielectric by means of its frequency dependent 
dielectric function e(k, u). This theory is not fully microscopic in the sense that the 
effects of the microscopic degrees of freedom inside the dielectric are embodied in the 
dielectric function but otherwise not explicitly described, and the Maxwell fields are 
subjected to the macroscopic boundary conditions at the surfaces of the dielectrics. 
Some insights on this theory and related approaches are given in the section 4.5.2. 
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Let us just quote here the specification of the Lifshitz formula to the short distance 
case corresponding to (151) 17 

f. f°° poo 2 

f(d) ~ — - / ds / du —————— (153) 

J 1 ; 16vr 2 rf 3 J J A A (iu)A B (iu)e s - 1 1 ; 

e A (uj) - 1 e B (Lu) - 1 

Higher order density contributions are now included in the dielectric functions. It is 
interesting to recover (151) by expanding the dielectric function to the lowest order 
in density 

e(uj) ~ 1 + A-npaiu), p — > (154) 
where a(u) is the electric susceptibility. Then A _1 (u;) ~ 2rrpa(uj) and 

/(d) ~ _hfA_PB_ f ^ ss 2 e ~s f duaA(iu)aB(iu) 
4<r Jo Jo 



fi-PAPB 

2d 3 



poo 

I duaA(iu)a B (iu) . (155) 
Jo 



This is precisely the expression (151) obtained in the elementary additive theory. 
Similarly taking the appropriate form of Lifshitz formula for long distance one can 
show that it reduces to (152) at low density as it should. 

It is remarkable that by the observation of the force (155) at the macroscopic level 
one can infer in principle by comparison to (151) the exact form of the microscopic 
van der Waals potential — CVr -6 for a single atom pair. 

4.4.2 Dielectric bodies at non zero temperature 

The Lifshitz theory can be worked out at finite temperature. In the high temperature- 
long distance regime characterized by a — > (21) one finds 

i r°° s 2 



I6ir(3d 3 J A A A B e s - 1 



87id 3 * — ' n 

n=l 

a e A + 1 a e B + l 

A A = 7, A B = 7 

e A -l e B -l 

where now e A = £a{w = 0), e B = e B (u = 0) are the static dielectric constants. 
In the perfect conductor limit of electrostatics e A ,e B — > oo, A^ = A B = 1 one 
recovers from (156) the classical Casimir effect f(d) ~ — 8 ^] 3 derived in section 3 

17 Here the dielectric functions e(k, lj) = e(ui) are assumed to be independent of k and magnetic 
polarization effects are not considered. 
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on a microscopic basis. We emphasize that Lifshitz has obtained the formula (156) 
by performing the large T asymptotics first, keeping the dielectric functions finite. 
The perfect conductor limit e — > oo is taken in a second step. Schwinger et al. [41] 
have proposed to take the limits in the reverse order, yielding the twice larger high- 
temperature force formula of Casimir (39). In view of our result of the microscopic 
analysis of section 3, we see now that the Lifshitz procedure is the correct one to 
recover the high-temperature regime for conductors. 

4.5 Theories 

All the theories are based on macroscopic Maxwell equations for the quantized elec- 
tromagnetic fields E, D, B, H in a dielectric medium. The system carries no net 
charge density and magnetic properties are ignored, so H = B. For oscillating 
fields E(x,t) = E(x,cj)e^*, B(x,i) = B(x, uj)e~ iuJt , D(x,t) = D(x, cj)e~ icJt , the 
equations for the frequency dependent amplitudes are 

VD = 0, VAE = i-B (157) 

c 

VB = 0, VAB = -i-D. (158) 

c 

One assumes that D is related to E by the linear and local relation 

D(x, oo) = e(x, w)E(x, u) . (159) 

The dielectric function e(x, a;) is piecewise constant in the bodies. From (157), (158) 
and (159) results the Helmoltz equation 

uj 2 

V 2 E(x, u) = -e(w)E(x, u), V • E(x, uj) = (160) 

cr 

valid in each dielectric domain characterized by the appropriate function e(u). The 
eigenmodes are determined by solving (160) under the usual boundary conditions 
for dielectrics: 

ei(u)E ijnorm (uj) = ej(u))~E jtnorm (u)) (161) 

for the normal component of the electric field at the interface between the domains 
i and j and all the other field components are continuous. 

4.5.1 Zero point energy method (Casimir's method) 

This approach is the simplest and follows Casimir's original idea that the force is 
generated by the modification of the zero point energy due to the presence of the 
dielectric bodies; it is presented in [4], sec. 7.2, and in [10]. Here e(u>) is assumed 
to be real which amounts to neglect the absorbtion of electromagnetic energy in 
the media. The new eigen frequencies are obtained from the zeros of the secular 
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equation, say R(u) = 0, determined by the Helmoltz eigenvalue problem (160), 
(161); they lie on the positive real axis of the complex cj-plane. They are generically 
noted {uj n } n and repeated as many times as required by their multiplicity. For 
two slabs at distance d, the eigenfrequencies uj n = uj n (d) depend on d and on the 
dielectric functions €a(uj), e#(u;). By the residu theorem the total energy can be 
written as 



E(d) 



J2^-huj n {d) 



1 

2ri 









-hw 




2 



d 

du 



\nR(u) 



(162) 



where C is the contour consisting of the imaginary axis of the complex w-plane closed 
by a large semicircle in the right half plane enclosing the eigenvalue. Then the force 
is defined by F(d) = —^E(d). Of course the non trivial calculation (which we do 
not present here) is that of R(uj) for given geometries and dielectric functions, as 
well as the discussion of various points concerning the convergence of integrals. For 
the slab geometry this leads to the same general formula as that of Lifshitz. Some 
special cases of it have been given above. 

When the temperature is different from zero, one introduces the free energy 
formula for each oscillating mode as in (22) (including the zero point energy) 



f(u) = -hw- /T 1 ln(l - e- pho ) = / 9~ 1 [2 sinh(/3^/2)] 



(163) 



and the total free energy <&(d) is given by (162) with [\hio\ replaced by f(ou). Notice 
that the calculation of the frequency spectrum is entirely classical and Planck's 
constant is introduced through (162) (apart from a dependence that can occur in 
an explicit expression for the dielectric functions e{uj)). 



4.5.2 Fluctuating fields (Lifshitz method) 

The Lifshitz theory [39] , [40] is formulated in the framework of stochastic electromag- 
netic fields. In the basic relation D(x, a;) = E(x, a;) + 47rP(x, uj), the polarization of 
matter P(x, uj) = P(x, u) + is supposed to be the sum of two contributions, 

a deterministic part P and a random part J£ which has zero average: = . The 
random part embodies all the quantum fluctuations effects arising from matter and 
fields. The relation (159) still holds for the averaged field 

D(x, u) = e(x, cu)E(x, to) (164) 

where the dielectric function is the same as before, so that for the fluctuating field 
(159) is replaced by 

D(x, u) = e(x, w)E(x, u) + K(x, uo) . (165) 

Then the macroscopic Maxwell equations in the form (157), (158) are supposed 
to hold only for the average fields E and B. Taking into account the random 
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polarization in (165), the actual random electric field generated from K obeys the 
equations with source 

V A E(x,w) = ^B(x,cj) 

V A B(x, u) = -i^D = -i^e(w)E(x, u) - i^K(x, w) . (166) 

The spirit is very much the same as in the Langevin theory for the Brownian motion 
where systematic effects of the medium (deterministic force and friction) are sepa- 
rated from the stochastic force (having zero average) due to microscopic collisions. 
As in the Brownian motion theory one has to define the correlations of the random 
force, due here to the zero point field fluctuations. They are taken of the form 



^(x^l^x'X) = 2hlm[e{u)]8 i> j8{u + u')8(x - x'). (167) 

Different Cartesian coordinates and frequency components of K are not correlated 
and the spatial correlations have zero range. This is equivalent to the white noise 
assumption for the stochastic force of Brownian motion based on the different time 
scales for microscopic collisions and macroscopic motion. Here the spatial scales are 
the intermolecular distances compared to the range of variation of the macroscopic 
Maxwell fields. The inhomogeneous Maxwell equations (166) are solved with the 
dielectric boundary conditions providing the fields E and B in terms of the random 
polarization K (together with the dielectric functions and the geometry of the bod- 
ies). For two slabs at distance d, the force per unit area on the first plate is given 
by the xx component of the Maxwell stress tensor T xx (x, c^)| x =o integrated over all 
frequencies and averaged on the random polarization according to (167) 



f(d) = - I duT xx (x,u)\ x=0 (168) 

where 



T =J- 

xx 4vr 



^ + ^-l(|E| 2 +]B| 2 ) 



(169) 



When the temperature is different from zero, (167) is replaced by the statement 
of the fluctuation-dissipation theorem (taking into account quantum and thermal 
fluctuations) 

K i {K,u)K j (x , ,u') = 2hTm[e(u)]6ij6(u + a/)£(x - x') coth (170) 
which reduces to (167) as T — > 0. This is the way h and T enter in the theory. 
4.5.3 Source theory (Schwinger's method) 

J. Schwinger, together with a number of physicists, were dissatisfied by derivations 
of the Casimir force relying on the mechanism of vacuum fluctuations. He proposed 
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to derive the results of the Lifshitz theory from a "more sound" basis where electro- 
magnetic vacuum fluctuations are not invoked [41]. The description (following [4], 
sees. 7.4, 7.5) starts now at a microscopic level by noting that a dielectric medium 
is characterized by a microscopic polarization density P(r, t) = ^2i^i(t)5(r — !•$), 
where Vi(t) are the induced dipole moments of atoms located at positions iv These 
dipoles act themselves as a source of a (microscopic) electric field obtained by solving 
the Maxwell equations yielding 

E?(r,t) = 8n J dr'dt'G^ir, t\r', t')P a (r', t') (171) 

where G is a Green function. The total electric field Eo + E s is the sum of a free 
quantum field E plus the field E s due to the dipole and the total energy 

U = -\j dr(P(r, t) ■ E(r, *)) = ~J dr((E+(r, f) • P(r, t) + P(r, t) ■ E"(r.t))) 

(172) 

is understood as the average energy of these dipoles in the field that they themselves 
produce. Here E + (r, t), E _ (r, t) denote the creation and annihilation parts of the 
field (see (137)). The average is taken on the vacuum of photons and on the state 
of the atomic variables 1S . The order of operator in (172) is in principle irrelevant 
since matter and field operators refer to independant degrees of freedom and thus 
have zero equal time commutators. The choice of order is however not innocent in 
subsequent calculations. Here the choice is normal order, placing the annihilators 
to the left of the creators. It is then evident that there is no contribution of the 
vacuum field in (172) since a^O >= 0, < 0|a£ A = 0. Then inserting the positive 
and negative frequency parts of (171) in (172) gives 

U = SnRe^J dr J dr' J dt'G + '^(r, t\r', t')(£"(r, t)P a (r', t')) j (173) 

with G + the positive frequency part of the Green function. 

To continue the calculation one has to know the atomic dipole correlation func- 
tion (P fl (r,t)P a (r f ,t')). For this one makes the simplifying assumption that atoms 
form perfect gases of uncorrected identical dipoles. To obtain the force one must 
further adapt the formula (173) to the two slab system and evaluate the Green 
function in this geometry. Within this scheme one eventually arrives at the Lifshitz 
main formula without any intervention of vacuum fluctuations. 

Other types of ordering can be used in (172), e. g. symmetrical ordering where 
the vector potential is not split in positive and negative frequency parts along the 
calculation. In that case the Lifshitz force can be seen to arise exclusively from 
vacuum fluctuations 19 . The final form of the Lifshitz force is of course indepen- 



18 In a stationnary state U will be time independent. 

19 As an example, the computation of the retarded van der Waals forces from vacuum fluctuations 
in section 4.3 involves the symmetrical order, see (142). 
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dent of the ordering convention, and can therefore receive different interpretations 
according to different choices and to the taste of the readers. 

5 Casimir effect in critical phenomena 
5.1 Introduction 

It is well known that at a critical point or in a phase with broken continuous symme- 
try, a statistical mechanical system in its bulk phase exhibits long range correlations 
(decaying algebraically rather than exponentially fast). This can be interpreted as 
generating a long range effective potential in the system. If this system is now 
confined to a slab of thickness d (e. g. a liquid film), then the development of 
these long range fluctuations as d — > oo may be the source of a force of Casimir 
type between the faces of the slab (film). More precisely the effect is reflected in 
the behaviour of the finite size corrections to the thermodynamical potential in the 
confined geometry. 

Let 4>L(d,T) be the thermodynamical potential for the slab of surface L v ~ x and 
thickness d in dimension v and 

<p(d,T) = lim ^0 L (rf,T) (174) 

L^QO Li 1 

the corresponding thermodynamical potential by unit surface. The finite size scaling 
analysis in the critical regime shows that is large d asymptotics has usually the form 

<p(d, T) = cVbuik(T) + yw (T) + + • • • (175) 

Here y?buik(^) is the potential density in thermodynamic limit and ^p SUT {(T) is the 
surface potential correction. The coefficient A(T) in the third term in (175) is called 
the Casimir amplitude. The Casimir force between the slab faces is defined as 

f(d) = -d d [^d,T)-d Vhulk (d,T)] = (v-l) ^p- + ... . (176) 

At a critical point or in phases with long-range correlations generated by the broken 
symmetry, the value of A(T) is expected to be non-zero and universal, depending 
only on the system and the boundary condition universality classes. Out of the 
critical regime the finite size corrections are expected to be exponentially small, 
thus the Casimir amplitude A(T) = 0. The tractable quantum models with phase 
transition where these ideas can be checked are not so numerous. A candidate, the 
free Bose gas that shows the phenomenon of Bose condensation, will serve as an 
illustration. For a thorough discussion of the subject, see the books [12], [13]. 
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5.2 The free Bose gas 

Here we report on a joint work [22] with V. Zagrebnov. We consider a free Bose gas 
in a slab of thickness d with faces of surface L 2 and periodic boundary conditions 
in all directions. As it is well known from the standard treatment of the free Bose 
gas, the finite volume grand canonical pressure $l(c£, T, //) at temperature T and 
chemical potential \i < is 



GMT,?) = i 5>[1 - exp(-/3(e(k) -/.))] 



;i?7) 



with the wave numbers and energy given by 



2™x , 



27m y 27m 2 

7 5 = 7 
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2m 



2im x 



2 7m, 



The potential per unit surface is then 



n x ,n y ,n z G Z 



27m. 



(178) 



¥>(d,T,/x) = lim -2* L (d,T,/x) 

L^OO Li 



= T2W0 1 dq £ in f 1 - exp (-^ - 2 " 2a2 G) 2 + ^) 

^ ' ' n=—oo L ^ ' 

(179) 

a 2 2 

where q = (q y ,q z ),Q = is a two dimensional wave vector, (3e(q) = and 
A = h\J [3/m the thermal wave length. Writing in radial coordinates J dq ■ ■ ■ = 
2n / °° dqq - ■ ■ and performing an integration by parts one can also write 



(p(d,T,fi) 



' n=—oo tJXJ 



(3e(q) + 2n^ {i)-0lM 
1 



- 1 



2%{3\ 

00 



exp 



+ 2 n*\* (f) 2 -/5/i 



(180) 



;isi) 



\ 2 2 

where we have introduced the dimensionless variable v = (3e(q) = and the 
function 

1 f°° 1 
^ U) = -2^J rf " exp[n2W-^]-l (182) 
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5.2.1 The normal phase 



The regime \x < caracterizes the normal phase (absence of Bose condesation). In 
this regime, the function ip(u) is infinitely differentiable so that we can represent the 
sum by the Euler-McLaurin formula. Since ip(u) = ip(—u) and using (18) one has 

oo oo 

^d- 1 n) = 2Y,*P(d- 1 n)+m 



= d 2 



POD 

/ dwp(u) 
Jo 



n=0 



^a» (0) + |l,(3» ( o) + ^» ( o) + ..: 



(183) 

Coming back to the definition (182) of ip(u), one sees that the first term of the large d 
expansion (183) equals -d PbuikiT, fj) where p hu i k (T, jj) = - lim^^ lim^^ ^^ L (d, T, jj) 
is nothing else than the bulk pressure, as it should. Clearly all odd derivatives of 
the even function ip(u) vanish at the origin, implying that finite size corrections to 
ip(d,T,ii) vanish faster than any inverse power of d, 

(p(d,T,fjL) = -dfr^T, ii) + O(j) for all r > 0, /i < (184) 

One therefore concludes that A(T, //) = 0, /x < in the normal phase (in the 
particular case of periodic boundary conditions, there are no surface terms). 



5.2.2 The condensed phase 

Bose condensation occurs in the free gas when the chemical potential is set equal to 
zero. Then the situation is very different : when /i = derivatives of ip(u) diverge 
at the origin and one needs a different method. We first expand for fi < the 
fraction in (180) (the Bose distribution) in power of the activity e 13 ^ and perform 
the ^-integral. This yields 

r=l n=— oo x ' 

The n-sum can be dealt with a version of the Jacobi identity stating 

oo 
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so that, setting a = * 2 r , we obtain the following exact representation of the slab 
potential 



(p(d,T,fi) = 
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(187) 
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since the first term yields the familiar low fugacity expansion of the bulk pressure. It 
can be checked that when p < the finite size correction term is (9(exp (—^—pd/2)), 
thus giving a precise exponentially small estimate of the remainder obtained from 
the Euler-McLaurin formula (184). 

For p = the r-series in (187) is convergent and we can write it in the form 
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1 oo 1 
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(188) 



The value of the integral is and when this is inserted in (187), one obtains the 

final result 

<p(d, T, p) = -d p bu i k (T, p) - + O (189) 

which yields a non zero Casimir amplitude A(T, p = 0) = — according to the 
general definition (175). 

One can impose different boundary conditions on the faces of the slab, e.g. 
Dirichlet, k = n ( n ~f~^ ; Q r Neumann, k = n = 0,1,2,... and find out by the 



same methods the value — |M for the Casimir coefficient in these cases (there are 
additionnal surface term contributions which are not present when periodic bound- 
ary conditions are used). 

Now few remarks and comments are in order. 

The grand potential of a free Fermion gas does not have a Casimir term for any 
value of the chemical potential. Indeed replacing the Bose by the Fermi distribution 
in (181) gives a corresponding ip(u) function (182) that is infinitely differentiable 
at u = for all /i. Since the function and all its derivatives vanish at u — the 
Euler-MacLaurin formula always yields corrections smaller than any inverse power 
of d. 

We note that the Casimir terms found in (189) are classical and universal, namely 
they do not depend on the Planck constant and the particle mass. In fact, for the 
free Bose-gas, it follows from a simple dimensional analysis that such a term is 
necessarily of the form C ksT /d 2 where C is a numerical constant. It is present at 
all positive temperature provided that the density p of the gas is higher than the 
critical density p c (P)- According to common wisdom, the Casimir force is due to 
Goldstone modes (i.e. low energy excitations) that will occur in the bulk limit when 
a continuous symmetry is spontaneously broken. In the grand canonical free Bose- 
gas it is explicitly seen that the excitations e(k) — p become gapless when p = 0. It 
is known that this generates long-range particle-particle correlations 



p( 2 )(r 1 ,r 2 )-p 2 ~p (T)|r 1 -r 2 | 



(190) 
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as |ri — — > oo in the condensed phase ( London- Placzek formula. [42]). Here 
po denotes the Bose-Einstein condensation density of the perfect Bose-gas. Casimir 
forces are usually attributed to such correlations in the critical regime. 

In the case of electromagnetic interactions, the Casimir term is always present 
as a result of the long range of the forces. In the standard calculation of the zero- 
temperature Casimir force between perfect conductors (section 2.2) the Casimir term 
appears because the third order derivative (19) occuring in the Euler-MacLaurin 
expansion does not vanish, contrary to the case at hand (183). This is due to 
the linear form fkv^ = Hc\k\ of the photon spectrum (non analytic at k = 0). 
Massive photons fko^ = c^J (ft|k|) 2 + (mc) 2 do not produce Casimir forces for the 
same reasons as for Fermions. 

We have studied in detail the classical electromagnetic Casimir effect in section 
3. In the spirit of the present analysis, one can say that such Coulomb systems are 
critical at all temperatures because the potential fluctuations are always of a long 
range, although charge correlations are themselves of a short range as a consequence 
of screening. 
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